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Abstract © 


The transient response of a dipole antenna is investigated using the singularity expansion method (SEM). The 
required natural modes of the dipole are obtained to good approximation using the finite difference time domain 
(FDTD) method. In order to generate a certain mode, a plane wave is used to excite the dipole antenna with both 
a frequency near to the natural dipole mode (resonance) frequency and with a spatial distribution approximating 
the mode dipole current. The method is applied to a dipole loaded with a resistor, thus the SEM solution is 
obtained without the need for the dipole impedance in a wide frequency band, The transient response of a dipole 
is studied for an incident nuclear electromagnetic pulse (NEMP) for dipoles of different lengths, and for an 
incident short damped sinusoidal plane wave, representing a high power microwave pulse, with different 
frequencies of the incident wave. 

Introduction. 

Investigation of the transient response of antennas is important in studying EMP effects, high power 

microwave weapons, localized waves, ultra-wideband communication and radar systems, EMP sensors ,..etc [1- 
6]. Analysis of the transient response of different types of antennas has been implemented-using time domain 
integral equations (IE), frequency domain IE together with Fourier transform, SEM, FDTD, in addition to high 
frequency techniques such as time domain physical optics, time domain geometric theory of diffraction and time 
domain puised plane waves [7-12]. 
The SEM has the advantage that once the antenna natural modes are obtained one can solve different excitation 
problems in a simple and fast way using a few numbers of natural modes [13]. Such transient problems include 
exciting the antenna with different waveforms (either as a transmitting antenna or a receiving one) and with 
different angles of incidence of a plane wave. The solution is particularly simple after the early time period 
when the incident wave has traversed the antenna. 

The natural modes complex frequencies correspond to the resonant modes of the antenna with damping 
factors accounting for antenna radiation during current oscillation. However, obtaining the natural modes is 
usually cumbersome, based on solving homogenous integral equations numerically and searching for the 
complex natural frequencies [14]. 

In the present work the basic natural modes of the dipole having minimum damping, which are required for 
solving transient problems [15,16], are obtained in a simple way by using the FDTD method. The modes are 
obtained by exciting the antenna with a transient plane wave containing a frequency spectrdm concentrated near 
the frequency of the required mode and with a spatia! field distribution similar to that of the required mode 
custent. After the transient excitation has died out, the excited antenna oscillates in the intended mode with an 
exponentially damped sinusoidal waveform, which is used to find the mode complex frequency. By such proper 
choice of the temporal and spatial distribution of the exciting waveform, each of the resulting modes becomes 
highly pure, i.e. least contaminated by other modes. This method can be used to simply obtain natural modes of 
dipoles with loads, which simplifies the solution of the transient problem without the need to have information 

_ about the short circuited modes of the antenna and the antenna impedance in a wide frequency band [15]. 

Once the antenna modes have been obtained, it is straight forward and simple to obtain the transient response 
of the antenna to different excitation waveforms and different incidence angles. Transient response of dipoles is 
then investigated for NEMP and short pulse waves resembling high power microwaves with good accuracy. 


The SEM method for dipole antennas 
The electric field integral equation for the current I(z) on a dipole antenna with exciting electric field E;, 
parallel to the dipole, is [16-18], 


LI = -sé,E, (1) 


Here L is an operator given by 
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where the dipole occupies 0<z<I, s is the laptace transform variable, c is the speed of light, &, is the free space 


permittivity and R is the distance between a current point at z' and an observation point at z on the dipole 
surface. 


For a dipole of length 'I' and radius 'r’ the antenna parameter Q2 is given by 
i 

Q=2In(-) (3) 
r 


For a thin dipole with large © the asymptotic natural modes currents I,(z) of the dipole, obtained from the 
solution of the veces version of equation (1), are sinusoidal, given by {16] 


1, (@)= sin) ; a (4) - 
The complex ine frequencies s, are given by 
3, = J Wr * oC, (5) 
=Inncl-o, 


where the damping factors 0, are proportional to 1 / () and J is the imaginary unit. For dipoles with QQ which 


is not large, the complex natural frequencies should be obtained numerically {14,15]. For an incident transient 
wave on a dipole, the transient excited current on the dipole can be obtained in SEM method mainly in terms of 
the natural modes of the dipole, with additional polynomial and entire functions of s, whose contributions are 
usually neglected for dipole antennas [16,19}. The Laplace transform of the excited current is given by [16, 19] 


I(z,s) =-€)s),(s—5,)' <El,,1, >“< E,1, >1, 6) 
where the summation is over the complex frequencies s, and their complex conjugates, and the operator 
L =dL/ds 
and the scalar product is defined by 
< fg >= {'S(@.sz)e (7) 

Thus 
—€)S5 ¢ 
I(z,s)= 9), ——1,(z (8) 
(z,s) DB ar pa) 
where c, is the mode coupling coefficient, given for the thin wire natural modes by 
i 
e, = JE.(,s),(2) @) 
0 
And B, is the mode normalization constant 
B,=<L'l,,1,>" 
= (ey (10) 
4c 


The Laplace transform of the excited current thus takes the form 


I(z,8) = oe Zo = 5 fin LE, 2" s)\dz ain) +000" (11) 


a 


The dipole current for an incident plane wave with damped sinusoidal time 


dependence 
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Consider an incident plane wave with damped. sinzsoidal time.dependence, piven.by 
E(t) =2E, cos(w, te (12) 


where wy is the radian frequency and o, is the damping constant of the incident wave. The Laplace transform 
of E(t) is. ee E,=1) 


1 {13} 
E(s) = ——+ 
§-a s- a 
Where © denotes the complex conjugate and 
a=Jw,-o, (14) 


The damped sinusoidal wave can be used to represent an high power microwave signal and can be reduced to 
the EMP signal. 


if the plane wave is obliquely incident on the dipole with an angle 6 with the z- axis, the wave front arrives at a 
point z on the dipole after a time delay 


zcos(@ 
eee) (15) 
¢c 
This time shift leads to multiplying the Laplace transform of the incident wave by e~” , thus 
1 | & Cun 
B(s)=(——_ + —e" = } (16) 
S-a s-a gieaa' £— @, 
The nore ae due to E, , thus equals 
sl_(z)<e "I > 1 
M(z, )=5 hi 2) en a sre ce (17) 


8 (s-s,) 4,74,0" S—a, 
To obtain oe time dependence of the current, the inverse Laplace transform is performed. By using partial 
fractions 
S, a, 
ae (18) 
(s-s = a,) ~(s- 5,)(s; - a;) "ee a,)(a, - s;) 
‘The required inverse Laplace transform for a general term in eqn. oD) is thus 


3 
r' e- ett) 


s 5 aed 
sohget ora; ,and the response starts at =T. Thus 
= (-T) 
L(z)a| 5, se ”1,> se alae 
I y= 5 int a of 20 
@t) in(®) fae a, Per ae (20) 
eet 
f (o)= [e* sin Ede (21) 
where T is a function of z given by eqn.(15), then 
Ay 1 AS i 
f,(0) = {| —————— |-| 2) (22) 
cos@ 
-bcos(@)/e+ A, =F - oli, 


where A, = J a . The coupling coefficient function in eqn.(21) is not a function of s (type 1) [22]. For an 
incident field of amplitude Ep in eqn.{12) rather than 2Ep , eqn.(20) reduces to 
21cE, et 

Cio Fie DS I(z) ¥ |= Fs) +4, a ie| (23) 


SM Se a;*aa" f a; if i 


It, nae 
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The exponential terms including s, represent natural oscillations of the dipole, whereas those including a; 
represent forced oscillations. 


Calculating the natural modes of the loaded dipole using FDTD Method 
“The FDTD method is used to find both the natural dipole modes with a resistive load and also to check the 
dipole current for an incident transient plane wave. The conventional Yee algorithm is used to obtain the 


scattered field in free space [9, 10]. Liao absorbing boundary conditions are used to enclose the computation 
domain. 


FDTD equations at a cell with a lumped resistor 

We consider a dipole antenna loaded with a resistor R. The FDTD updating equations at 
the cell with a lumped resistor R, oriented in the z- direction, can be obtained by relating the current density j, 
flowing in R to the total parallel electric field in the cell using Ampere's curl equation for the total field (20) 


i) ae 
Vx H, =€ ea +j (24) 
‘The total field is the sum of the incident field E, and the scattered field E, at a cell (i, j, k), for example 
E, (i, j,k) = £,(, j,k) + £,G,j,4) (25) 
In order to obtain the equation for the scattered field one considers the equation of the incident field. 
OE 
VxH, =e Sai _ (26) 
By eliminating the incident field from equation (24) one obtains for the scattered field 
OE, 1 
—=(VxH,-j,)— (27) 
ar ( oe 


whose finite difference form reads for the component E, at time step n+1 (with subscript 's' omitted), 
ALD Eee ; 
At 
eal all 1 
ards yi J,k)-H aG- Lik) HOG i K)- HG I-L’ TL, 274) 
é Ax Ay AxAy 


(28) 
where j,=I/A and I is the total current evaluated at time step (n+0.5), similar to the magnetic field, Ax, Ay, Az, At 
are the spatial and temporal steps (spatial steps are equal) and A is the cross sectional area 


A=Ax, Ay 
By eens the current in terms of the voltage 
PS=v/R= Az E™?, /R (29) 


The value of the electric field at time step n+.5 may be obtained as the average of the values at time steps n and 
n+ at which the electric field is calculated, thus 


EG 5,K)- Ee I) | 
AE 
HH"), j,k)- Hy G-LiW) AG, ,-H XG j-1.) 
1 Ax Ay 
e|_ af GAB GIK , EN GIDTE ih »| (30) 
2RAxAy 2RAxAy 
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E™ , j,k) 14—\—_1= E."(i, j,k} 1 -——_ | + —- 4, 06,3,0-A,0-L), 
te er | "i | “5 ft 6, 5,4)-H,G-L 4,4) 
“At AzAt 
-~—[H,(,j,0-4,(6j-1L0)|- JEG 3,0" + E04)" 
ay 4B) ij Ram BiH) +E,(,j,b)" | 
(1) 


Excitation of the natural modes of the dipole 


The natural damped sinusoidal modes of the dipole take the form 

I,(t} = i, cos(co, t) exp(-c, t) (32) 

Where on and g, are are represented in eqn.(5). The main modes required for the SEM solution are the least 
damped modes [15]. Other natural modes exist with higher damping factors, but such modes are not relevant in 
the SEM. 

To excite a certain mode on the dipole, a plane wave is used to excite the dipole with approximate frequency 
and spatial distribution of field along the dipole as the mode current required to be excited 


E, = E, sin(1Z > sin(@,t)sin(w,t/8) (33) 


Where n is the mode number. In order to limit the period of the incident plane wave, a few cycles of the natural 
frequency o, are used, modulated by a half sine wave of a lower frequency (1/8 in eqn.33), Fig.(ta), in order to 
make the field continuous with time, thus the spectrum of this wave becomes concentrated around m,, so the 
other natural modes are only negligibly excited. After the termination of this excitation, the current oscillates on 

~-- the -dipole mainly: with the intended natural mode. The complex frequency of the mode is calculated by 
considering the successive positive and negative maxima P; and Pj,, in the damped sinusoidal current waveform 
and calculating the damping factor (0 ,) and frequency f, of the mode as 


1 
= — 34 
tn = Fal 
in P, -In Py Re 
> mat 


Where m is the number of time steps At between such two successive peaks. 


Results and discussion 


The natural modes of the dipole 
The short circuited dipole 
Fig(1b) shows the variation of the excited current, at the center of a short circuited dipole, with time, where the 
frequency of the incident plane wave is nearly equal to the frequency of the first natural mode, eqn.(5), n=1, and 
the direction of the incident plane wave is normal to the dipole axis. The forced excitation due to the incident 
wave terminates after 4 cycle then the damped oscillation of the first natural mode begins. The dipole used is a 
strip dipole of rectangular cross section with thickness to width ratio of 1/3 which is equivalent to a cylindrical 
dipole of radius to width ratio of 0.38 [21]. The dipole is 39 cells long and 3 cells wide, thus its equivalent 
radius is 1.14 cells (thus the corresponding Q = 7). The actual used frequency is 6 GHz, with Ax=0.000625m 
’ and At=1.2*10 sec. This time step equals the cell length divided by c*sqrt(3) to satisfy the stability condition 
[9,10]. A load of 1 Ohm is used in place of the short circuit to avoid overflow in the computations, eqn.(31). 
Fig (2a) show the current distribution of the first mode on the dipole, which is seen to be nearly sinusoidal. Fig 
(2b) shows also the current distribution of the third mode which is also nearly sinusoidal, showing nearly no 
contamination by the first or fifth modes (which have even current distributions w.r.t. the dipole center), The 
even numbered modes contain an even number of half sine waves and have zero current at the center of the 
dipole. 
Table (1) shows the computed damping factors and netural frequencies (rad/sec) of the first 5 modes, 
normalized w.r.t. twice the nominal first resonance frequency of the dipole (c/l). Comparison is made with the 
results of [15]. The comparison shows reasonable agreement, particularly for w,, with the deviations expected to 
be due to slight contamination by the other natural modes and thé small difference of the dipole parameters 0 of 
the 2 compared cases. 
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Dipole with 75 Ohm load 

Fig(3a) shows the current distribution of the first mode for the dipole loaded with 75 Ohm. Fig (3b) shows the 

ae distribution of the third mode. A slight dip can be seen in the current at the dipole center due to the 75 
load. 

Table (2) shows the calculated damping factors and natural frequencies. The damping factors are higher for the 

odd numbered modes for the 75 Ohm case than for the 1 Ohm case (table 1) due to mode attenuation by the 

resistive load. On the other hand, for the even numbered modes the damping factors are nearly the same since 


the current is zero at the dipole center for these modes, thus the damping is not affected by the presence of the 
load. 


Table (2) The modes normalized damping factors and radian frequencies a for dipole loaded with 75 Ohm 
damping factor, o(V/c) 


Sa eee ba eee 
Saar eee | eee 
3 as a ee, 
Sree LE Lien FE i eee 
LSS fe eae ae 


Dipole response to an incident NEMP 
An NEMP can be approximated by the double exponential waveform [15], [2] 


E=E,(e -e*) (36) 
where @ = 1.6*10% sec’, B =2.6"10* sec"! 


B Corresponds approximately to the inverse of the rise time constant, and @ corresponds approximately to the 
inverse of the decay time constant, and E, equal aapproximately 5*10° V/m for NEMP. The laplace transform of 
Eis 

Gh Roe 
esta s+B (37) 

Each of the two exponential factors can be considered as a special case of a damped sinusoid with zero 
frequency, eqns. (12, 13). 

Fig.(4a) shows the computed transient current in a short circuited (1 Ohm) dipole of length 2 meter (75MHz 
resonance frequency) due to an NEMP normally incident on the dipole axis. In this case the even natural modes 
are not excited. Three odd modes are used in the calculations, which are enough to get good accuracy at 
relatively late time. The main excited mode is the first mode. Figure (4b) shows the computed response using an 
equation of Ref. 15 corresponding to ¢qn.(23), with the exact modes resonance frequencies and damping 
coefficient taken from [15] (table1). Fig.(4c) shows the computations of the excited current using FDTD method 
for an incident NEMP. The comparison shows good agreement of the 3 results. 

Loaded dipole 

The more practical case of a dipole loaded with 75 ohm is then considered. Fig (5a) shows the transient current 
for dipoles of different lengths resonating at different frequencies. It can be seen that the dipole current 
decreases as its resonance frequency increases. At high dipole resonance frequencies, Fig.(Se,d), the weak 
forced excitation by the double exponential NEMP waveform begins to dominate the current natural 
oscillations. The decrease of the current may be attributed to both the decrease of the effective dipole length (by 
100-times from 10MHz to 1GHz) and the decay of EMP spectrum, Fig.(6), eqn.(37), by approximately the same 
order of magnitude. The damping of the current due to the load can be seen by comparing Fig.(4a) and Fig.(5b). 
By using the loaded dipole natura! modes, one does not need to have both the short circuited dipole modes and 
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the dipole impedance data in a wide frequency range and does not need to use numerical inversion of Laplace 
transform as done in Ref.15. 


Incident damped sinusoidal! plane wave 
High power microwave weapons use few cycles of microwave frequency in order to penetrate devices operating 
at microwave frequencies. Here we use a damped sinusoidal incident plane wave to investigate such conditions. 
In order to investigate dipole response to an incident plane wave of short period (wide spectrum ), we use a large 
damping constant o,=1.9*10'°, which equals the inverse of the damping time constant of the sinusoidal wave 
envelope. If the frequency of the sinusoid is low, the high damping constant makes the waveform to be nearly an 
impulse. 
A dipole of length 2.5 cm (nominally resonating near 6GHz) is- considered with incident plane waves of 
different frequencies (f§,=1GHz, $GHz, 10GHz, 15GHz), normally incident on the dipole axis. Fig.(7) shows the 
dipole responses. In Figs.(7a,b), with {,=1 and 5 GHz, the first natural mode of the dipole is excited. However, 
in Fig.(7b), the second and third peaks are higher than in Fig.(7a) since the the frequency of the incident wave is 
near the dipole resonance frequency. In Figs.(7c,d), (f,=10GHz, 15GHz), contributions of the third natural mode 
are noticed, particularly in Fig.{(7d). It can be seen that the peak amplitude of the excited current varies slightly 
with changing the frequency of the incident wave. This can be attributed to the wide spectrum of the short 
period incident wave shown in Fig.(8), eqn.13. The peak excited current in a dipole is thus expected. to be 
“ maitily proportional fo its length when a short high power microwave pulse (with wide spectrum) is incident on 
the dipole, 
The response is computed for t >] cos@ /2c, when the wave arrives at the dipole center. The peak of the current is 
less than that of Fig.7b and the second natural mode (with odd spatial distribution w.r.t. the dipofe center) also 
contributes to the current (not shown at dipole center). Both of the effects are due to the inclined angle of 
incidence of the wave. ay 
Conclusion 
The computationally efficient SEM method is used to investigate the transient response of dipole antennas using 
the natural complex frequencies of the dipole. A FDTD approach is used to simply obtain the approximate 
complex natural modes of the dipole. The modes are obtained for a dipole loaded with a resistor, thus the 
solution of the transient problem does not require the dipole impedance to be known in a wide frequency range. 
The obtained natural modes are used to find the response of dipoles of different lengths to the NEMP. The 
excited current decays sharply as the dipole resonance frequency increases due to the decaying spectrum of the 
NEMP with frequency and the shorter effective length of the dipole. The response of a dipole to a highly 
damped sinusoidal plane wave (resembling high power microwave pulse) is obtained. The peak excited current 
in the dipole changes slightly with the main frequency of the incident wave due to its wide spectrum caused by 
its sharp damping. 
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Fig(1-a)The incident waveform 
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Fig.(i-b) The current at the center of a short circuited dipole versus 
time when excited by a plane wave with the waveform of Fig.{1a) 
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Fig (2-b) The current distribution of the third mode 
The current distribution of the first mode on the 
Short circuited dipole (the dipole occupies the cells from 
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Fig.(3-a) The current distribution of the first 
mode on the dipole loaded by 75 Ohm 
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Fig.(3b) The current distribution of the third mode on the dipole loaded by 75 ohm 
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Fig, (4) Transient current (Amp) induced at 

the center of a short circuited dipole of length 

2 m (nominal resonance frequency of 75MHz) 
with an incident NEMP (Ep=1000, At=10" sec.) 
a) Computed (1 Ohm load) 


Fig.{4 -b) Results based on Ref. 15, 0 Ohm 
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Fig.(4-c) FDTD results, 1 Ohm load 
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Fig.(5-a) Computed transient current (Amp) at the 
center of a loaded dipole with 75 Ohm due to NEMP 
(Eo=1000V/m, At=10''sec.)a) dipole of nominal 
resonance frequency of 1OMHz 
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Fig.(5-b) Dipole of nominal resonance 
Frequency of 75MHz 
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Fig.(5-c) Dipole of nominal resonance frequency of 200MHz 
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Fig.(Se) Dipole of nominal resonance frequency of 2GHz 
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Fig (6) NEMP spectrum (log-log scale) 
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Fig.(7) Computed transient current (Amp) at the center of a dipole of length 2.5cm loaded 
with 75 Ohm due to a damped sinusoidal wave (Eo=1000V/m, 

At =1.2*10°?sec.) 

a) §,=1GHz 
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Fig.(7b) f;=5GHz 
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Fig(7-c) 
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Fig.(7d) §=15GHz 
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Fig.(8) Damped sinusoid spectrum, linear scale 
(f =5GHz, 0, =1.9*10!'°) 
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